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FORMULAS FOR LONG PERIOD RADIATION PRESSURE, 
LUNAR AND SOLAR GRAVITATIONAL EFFECTS ON THE 
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By David Fisher 
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Greenbelt, Maryland 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


SUMMARY 


Long period formulations are given for satellites in 
which the perturbations due to solar radiation pressure, 
lunar and solar gravitational forces are of first order in 
magnitude. This occurs for the classes of scientific satel- 
lites whose ratio of reflecting surface to mass are .07 cm 2 
per gram or larger, and for satellites whose semi-major 
axes are between 12,000 and 55,000 kilometers such as 
communications satellites and satellites having eccentric- 
ities as high as 0.85. 

The method of canonical transformations is applied to 
obtain analytic formulations which may be added to existing 
long period perturbation formulas given by Brouwer or 
Kozai for the earth f s oblateness. 

A useful application of these results would be to incor- 
porate them into existing programs to improve the accuracy 
of orbit determination. 
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FORMULAS FOR LONG PERIOD RADIATION PRESSURE, 
LUNAR AND SOLAR GRAVITATIONAL EFFECTS ON THE 
MOTION OF ARTIFICIAL SATELLITES 

by 

David Fisher 

Goddard Space Flight Center 


INTRODUCTION 

As a result of examining the data derived from satellite observations it has been established 
that the sun and the moon perturb artificial satellites to a significant extent even when the satel- 
lites are relatively close to the earth. In Reference 1, for example, it is shown that the order of 
magnitude of the disturbing functions for satellites whose major axes are only 12,000 kilometers 
are of the same order of magnitude as those of the disturbing functions due to the earth’s oblateness. 

An unexpected effect on the orbits of satellites was that due to radiation pressure. In Refer- 
ence 1 it is shown that the perturbations due to radiation pressure were significant on a satellite 
of the Telstar n type with surface area to mass ratio of 0.075 cm 2 per gram. A ratio of this 
magnitude is often obtained in scientific satellites. Balloon type satellites such as Echo 1 have 
had ratios as great as 125 cm 2 per gram. 

An analysis of several distant satellites given in Reference 5 indicates that the expansions 
for the solar and lunar disturbing forces are valid at least for satellites whose major axes are 
55,000 kilometers. This covers the range of the U. S. and Soviet communications satellites and 
includes satellites of eccentricity as high as 0.85. In this report the long period perturbations are 
developed using the same disturbing forces as those in References 1 and 5 and consequently the re- 
sults should have the same range of validity both for solar and lunar gravitational forces and solar 
radiation pressure. 

For the accurate prediction of the orbits of satellites these solar and lunar perturbing forces 
should be taken into account. It is the purpose of this report to provide analytic perturbations in the 
elements of satellites. These may be readily added to the perturbations in the elements of satellites 
computed by analytic theories which take into account only the earth’s oblateness such as those 
given by Brouwer or Kozai, Reference 3. 

In this report the method of canonical transformations, Reference 2, is adopted in order to derive 
the perturbations in the required form. 
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GENERAL FORM OF THE DISTURBING FUNCTIONS 


A useful tool In finding perturbations by the method of canonical variables is to obtain a de- 
termining function s, which in turn is derived from a disturbing function R. The function R, Ref- 
erence 1, is a Fourier series of the form. 


R - R s + Q(a) 2 ^ A* cos a . . (1) 

i 

Q(a) is different for each type of perturbing force and is given in Table 1. All symbols are defined 
in the appendix. 

Table 1 


Q Functions. 


Type c r Perturbation 

Q(a) 

9 

3 _ A 

Long Period Solar Radiation Pressure 

- 2 F m a 

Short Period Solar Gravitational 

16 n e m e' a ’ 

Long Period Solar Gravitational 

3 n e%' 

16 (l ~ e 0 J ) 3/2 

Long Period Lunar Gravitational 

3 n e m e f 

(independent of moon's node) 

16 (l-ejy't* 

Long Period Lunar Gravitational 

3 n e 2 m e j 

(dependent on moon's node) 

16 (i-e e 2)3/2 a 


R s and A i are functions of L, g, h, e, and i ec , where L, G, and H are Delaunay variables, e and 
i ec are the obliquity of the ecliptic and the angle of the moon’s orbit with the ecliptic respectively. 

L, e, and i ec are considered as constants. The arguments a. are defined by 

a. = j. a + k.n + . 

j j, kj, and -C. are integers which may take on values of 0, ±1, ±2. 

t« is the argument of perigee, n the right ascension of the node, and \ 0 the longitude of the sun. 
The functions A k and the arguments a t for the Radiation Pressure, Short Period Solar Gravitational, 
Long Period Solar Gravitational and Long Period Lunar Gravitational effects of the moon are given 
in Tables 2, 3, and 4. 
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Table 2 

Radiation Pressure Functions 



Table 3 

Solar Gravitational Fourier Coefficients. 
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When the disturbing function R Is also dependent on the node of the moon n tc , r takes the form 

R - QOO^Vi^Vl . (2) 

01 

The functions A m and the arguments a m . are given in Table 4. 

The arguments a B _, are of the form 


a = j . a; + k, 0 + q . H 

m, 1 M i i e c 

The quantities j , and k i are integers which take on values of 0, ±1, or ±2 as in Equation 1. 
The quantities q m , take the values ±1, or ±2. 

In selecting terms for the solar and lunar gravitational functions the following two criteria 
were adopted 

(i) Terms only to order sin 3 e were kept 

(li) Terms involving the eccentricity of the sun or moon were not kept 

THE RADIATION PRESSURE FUNCTION 

The disturbing function for Radiation Pressure is given by 

R - F^rcosS (3) 

where 

F a const = -4.63 x 10‘ 5 dyne/cm 2 
A => presentation area of satellite 
m = mass of satellite 
r = radius of satellite = a(l - e cosE) 

s = angle between the radius to the satellite (r) and the radius to the sun (r ' ) 

Thus 


cos S 


xx* + yy * + it * 
rr 1 
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where by Reference 4 we have 


cos S = cos 2 2 sin 2 2 cos (f + ^ + fi + + cos 2 cos 2 ■j cos (f + w + ft - 

+ sin 2 ^ cos 2 "2 cos (f + v - ft + \ fe ) + sin 2 "J sin 2 y cos (f + ai - ft - 


- ■J sin i sin e cos (f + a + \ 0 ) + j sin i sin e cos (f + u - X 0 ) . (4) 

The long period part of R denoted by R p is found by integration with respect to the mean anomaly 
i t thus 


If we recall that 




iSd-t 


d-t = (l-ecosE)dE 
r cos f - a(cos E - e) 

r sin f - a /l - e 2 sin E 


we can readily derive the useful result that 


(5) 


1 

277 



(f + a.) d-t 


3 

~2 ae cos a ( 


( 6 ) 


where a 4 is an arbitrary argument independent of f . 

Applying this formula to the expression for R p we find that R p can be expressed by the Fourier 
series 


R 


a 




cos a k 


where the functions A ; and arguments a ( are given in Table 2. 


(7) 
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SHORT PERIOD SOLAR GRAVITATION DISTURBING FUNCTION 


The short period solar disturbing function is so called because it contains the argument 
the true longitude of die sun. However, no short period argument of the satellite itself appears. 

If we apply die criterion of keeping terms only 14 ) to order sin 3 « we find from References 1 
or 4 that f 


R e = ll Vi 


( 8 ) 


Both and a are given in Table 3. The argument a* has the form 


a i - 2K o + ii " + k i n 


j i = 0, ± 2 , k i = 0,-1, -2 


An important property of A Q \ is that it can be separated into the product of three functions. 


Thus 


A 0i " A 0i (e)A i (e)A i (i) . (9) 

The independent variables €, e, and i are the obliquity, eccentricity, and inclination respectively. 
The independent variable for the inclination i is of course not to be confused with the subscript i 
used as a summation index. 


THE SOLAR AND LUNAR SECULAR AND LONG PERIOD DISTURBING FUNCTIONS 


The solar and lunar disturbing function, representing the gravitational effect of distant bodies 
on the earth satellite have of course the same general form. Figure 1 shows a skef< h of the earth, 
moon, sun configuration in order to present 


some of the elements appearing in the dis- 
turbing functions. 

We start by writing a representation of 
the disturbing function of the moon, R e , as ob- 
tained from Reference 1 or 4. The disturbing 
function for the sun is readily obtained from 
R c . Then R c finally is separated into two parts 
in order to obtain the long period lunar per- 
turbations analytically. From Reference 1 or 


T (MOON'S ORBIT) 



Figure 1 
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s 


4 we have 


1 n e m e 


t 3 . 2 

4 ft - e J)3/2 

U"I Sln l c ) 

^1 - j i ) 


^ n ffi f x 

+ 16 JjZ e l) 3 /l ^5 e J cos 4 2 sin* i e cos [Xi'-to-tt) 

- 5 e 2 sin i cos 2 sin 2i c cos (f! e - 2 <v ~ ft) + 5 e 2 sin 2 i ( 1 - -j sin 2 i c ) cos 2" 

+ (l + -5 e 2 ) sin 2 i sin 2 i e cos ( 2 Cl c ~ 2 ft) + (l + e 2 j sin 2 i sin 2 i ffi cos (ft c - Cl ) 

+ 5 e 2 sin 4 -j sin 2 i c cos (2ft c + 2 co~ 2ft) + 5 e 2 sin i sin 2 sin 2i c cos (ft e + ~ ft)^ (10) 

This expression gives both the secular and trigonometric terms of the lunar disturbing 
function. 


THE SOLAR SECULAR AND LONG PERIOD DISTURBING FUNCTION 

From Figure 1 it is seen that by sliding the moon's trajectory T to a position coinciding with 
the ecliptic, the configuration of the sun will be obtained. Consequently the long period disturbing 
function for the sun is obtained by replacing i c with e, setting ft e equal to zero, and replacing the 

parameters of the moon by that of the sun in the expression for R c , Equation 10. When this is done 

the expression for R 0 is given by 

_ _ 1 n e m 0 a * f. 3 . , \( 3 . ,.\ 3 n o< a2 V. 

R e ’ 4 (! V 1 - 2 sin £ M 1 ’ 1 Slnl 1 ) + 16 (i_ “ 2 )S/ 2 L* A oi cos a i • (11) 

The first term of this disturbing function is a function of the non angular elements only. From 
this term die contribution of the sun to the secular motion of the satellite is derived below. The 
remaining terms of the expansion contain the arguments a. which in turn are linear combinations 
of u andn. 

The functions A 0i and the arguments a. are given in Table 3. 

The functions A Sl can in turn be written as the product 

A ®i = A ©» ( £ ) Aj («) Aj (i) . (12) 
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N?rfSBagnro: 




The functions A i (e) and A i ( i ) involving the satellite elements are the same as those in the 
short period solar disturbing function. 

In addition the arguments a. of the long period disturbing function Contain the same linear 
combinations of u and fi as do the corresponding arguments a,' of the short period disturbing 
function. 

Thus 


a ; - jjCd + kjH 

where j ; =0, ±2, k s = 0, -1, -2 exactly as in Equation 8. 

THE LUNAR LONG PERIOD AND SECULAR DISTURBING FUNCTIONS 

The form for R e given in Equation 10 is not suitable for ready analysis since the elements of 
the moon i e and n c are periodic functions in the variable n ec . However by referring to Figure 1 
we obtain the necessary relations between i c , and n ec to facilitate analytic formulations. Thus 
from spherical trigonometry it follows that 


cos i e = cos e cos i ec - sin e sin i ec cos 


sin i. sin fi- = sin i sin H 

C C ec ec 


(13) 


In the disturbing function R c we find terms of the following three types 


( 1 - sin 2 i c ) cos r 


(«e + d 


and sin 2i c cos (n c + r ) , where r is a linear combination of cu and n. 

If we now apply formulas (13) we find after some algebraic manipulations that 

(l - | sin 2 i c ) cos r = ( 1 - sin 2 i ec ) (l “ "J sin 2 e ) - j sin 2e sin 2i ec cos (r + H ec ) 


sin 2 e sin 2 i ec cos (r + 2H ec ) - -g- sin 2e sin 2i ec cos (r - O ec ) 


3 . 


+ g- sin J € sin J i <e cos (r- 2fi ee ) , 
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(2Vt) = (l- 


3 2 
2 


I sin 2 e + sin 2i sin e cos" 


cos 


( T + n .c) 


+ sin 2 i 


cos 4 — 


PftS ( T + 50 




+ 


sin 2 i ec sin 4 \ cos ( r ~2fie C ) 


sin 2i c cos (O c 


+ r ] - (l - -j sin 2 i ec ) sin 2e + sin 2i ec cos 2 j (- 1 + 2 cos €) cos (r + 0 ec ) 


- sin 2 i sin e cos 2 cos (t + 2ft*- ) “ sin 2i 


/ 1 + O rne r \ 


+ 


sin 2 i ec sin e sin 2 % cos (r - 2ft ec ) • 


(14) 


The next step is to apply Equations 14 to the expression for 1^, Equation 10, given above. 
Where this is done the expansion for R c , keeping terms only to order sin 3 e , takes the form 


4(1 


e Z ) 3/i 


(. 3 . , \l 

3 . 2 . \/ 

' 3 . 3 .\ 

U"7 sin e M 

1 - T sin 2 i e J( 

1 - ~2 sin * i J 


_3_ "£ a 2 y 

+ 16 (l-e s f>Z-^ 


cos a . 


n c m C a2 

16 (l - e e ) 3/2 



m 


(15) 


The first term on the right-hand side of Equation 15 is a function of the non angular variables 
only. Consequently from this term the contribution of the moon to the mean motion of the satellite 
is derived. 

The second term is a Fourier series whose arguments a. are independent of the node of the 
moon Q ec and are consequently identical with the arguments of the long period solar gravitational 
function. The arguments a. and the function A ci are given in Table 4. It is interesting to note that 
the relationship 




(l-|sin2 



(16) 


holds. 
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The last Fourier series in Equation 15 contains terms with the argument a m . where 

a = j.aJ+k.O+q . f) . 

m, i J i i ec 

We also have the separability relation 

A = A . (e, i ) A • (e)A. (i) . 

m,i m ,» \ ec I i V / i V > 


(17) 


We have kept the subscript i to preserve the correspondence between the Fourier series of 
the solar and lunar gravitational expansions of the disturbing functions. 

The functions A m * and the arguments a m . are given in Table 4. 


THE MEAN MOTIONS 


If R s is the secular part of the disturbing function of the sun or moon then the contribution to 
the mean motions by the sun or the moon are given by the formulas 


l = 


* R s 

TT ’ 


TS ’ 


n = 


*k s 

nr 


Performing the differentiations indicated above the combining with the secular motion due to 
the earth’s oblateness given in Reference 1, we have 


l = 


f 3 ^J 2 a e 2 /L 3 \/ H 2 \ 3 ^J 2 ia ' 4 /L 7 \[ G G 2 

"Y " ^ L (g 3 )! 1 " 3 G 2 ) + C \G^/ j_ _ 15+ 16 T + 25 


/ G 

G 2 \ 

H 2 

/ G 

G 2 N 

i H 4 1 

(30 - 96 r - 

90 u) 

+ 

(105 + 144 f + 25 

L 2 ) 



45 3 4 a * e2 ^L 7 \/ H 2 H*\[ 

- m — r — / 1 3 ' 30 + 35 ^)j 


(i 


(l- y sin 2 i ec ) 


( l_e c) 


3/2 
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The symbols employed here are defined in the Appendix. It is shown below that differential 
coefficients derived from the mean motions, a^/dL, aiVaL, G(<j<i/dG) , G(dO/dG), G (a^/ait) , and 
G (dH/dH) are useful in forming analytic perturbations. Performing the indicated differentiations 
we now give the formulas for the above differential coefficients. 

Formulas for ac/av and ao/at 


5T 


1 f 9 J, •.» L 4 f , . 3 J 2 Ja e 4 L*f/ G G 2 \ 

IT & U- 5 c °s i) + 128 ~ ~T^ [(l 05-96 r - 125 — ) 

+ 270 + 768 § + 630 jjj cos 2 i - (ll 55 + 1440 £ + 225 jj) cos 4 ij 

U. 

+ m ^ [( 63 - 45 $) + (- 810 + 630 T 2 ) cos2 1 + ( 115S ' 945 ^) cos4 *]} 


3 ( 3 . . \ n e m o / 3 . , . \ n c m c 

- vi- 2 sini£ ) ^ + U" i sinJ *.J ^ . e 7 ji 


( 1 _ e 2 - io cos 2 i) 


12 


1 [9 j 2 a . 2 L 4 . 3 


et _ tt n 


a 4 G 4 


— COS 1 + 


32 


Jj 2 a e 4 L 8 17 G G 2 \ 

— ^8 [(l5-48 r -45pJco S 


+ (l05 + 144 £ 4 25 g) cos 8 i] + ^ e 2 (3 - 7 cos 2 i) cos ij 


■ (l“ 4 si " 2 £ ) 

Formulas for o> Q , fi 0 ,i H , 


W. 3 . , . \ n e m c 


(7 + 3e J ) 


cos 1 


[3 J 2 a e 5 L 4 / , . 3 J 2 2 a e 4 L 8 \t 
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\~ 168 + 54 — J + 36[75- 28 — j cos 2 i + 30^- 154+63 g ) cos* i 


128 a4 


a 4 G 8 



— 
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THE DETERMINING FUNCTION 

In this report the method of canonical transformations is applied to obtain perturbations in 
the elements. In the notation of Reference 5 we wish to determine expressions for the unknown 
variables L', G', H', V, g', h' in terms of the variables L " , G", H", l" , g", h" which are known 
(usually derived from experimental data). L”, G", H" are constants while l", g” ,h" are linear 
functions of time. The disturbing functions described above should be considered to be expressed 
in terms of the single primed elements. 

By extending the method given in References 2 and 5 we may write the determining function s 
as 


A 4 sin a. 

S = L " + G" g' + r! h' + Q(a) 7 7 


Hie functions Q( a), A. , a, are functions of the double primed elements while the a i are functions 
of the single primed elements. The expressions for 0(a), A., and a. come from the expressions 
for the disturbing functions and are given in Tables 1-5. The quantities a. are linear combinations 
of the secular motions which have been discussed above. 

THE PERTURBATIONS IN THE NON ANGULAR ELEMENTS 
From the definition of the determining function we find 


SL 


dS T 

Tv ' 
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To simplify the nomenclature we replace g by u, h by n, and drop the primed and double primed 
notation for the remainder of the analysis. 

From the representation of the determining function given above we find 

SL = 0 

SG = Q(a) j j r - coi a, 

SH = 0(a) ^ ' kj cos a. . 

The relation between the Keplerian elements and the Delaunay variables is given by the formulas 


Consequently we find 


G = yiua (l - e J ) 
H = G cos i . 


he = 


- S G 

eC^ sc 


Si = cTTHT f 008 ‘ SG " SH ] 


THE PERTURBATIONS IN THE ANGULAR ELEMENTS 
From the theory of canonical transformations we have 


hi - 


«9S t 

1l 


hcu - - 


ha 


dSj 

Jg 


s s T 
1h 


Where hi = V -i * 9 bu> - <*' sn = n # -n* and 


s t = 


16 



In order to find die perturbation In the mean anomaly It Is convenient to replace Q(a) by Its 
factors 


0(a) = Q' f(a) 


where Q' Is a constant. 
The as T /(?L Is then 


£fr 

<9L 


.. V ll lJ ( f (“) *i) 

Q L 


(f(a)A t ) 


'll 


*in cl : 


It then follows that 


U 



1 *»\ 

I-TTj «»“i 


where 


i 

u i • f(a) Aj <?L 


i»L 


, do; dfi 
Ji <JL + k i <»L 


where for example In the case of radiation pressure 


f(a)A i x factor independent of L 


for all i so that 


1 + e» _ He a 
e*L na J e 1 


for all i . 

In a similar manner It is found that 


hCA) 


. Q(«) 

' " G 




sin a 4 
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G dA i 
K i " A, <9G 


fi, 

a iO ' G .9G 


G 


du ah\ 
dC + k i dO) 


t.n = 




iin a 4 


G l9A i 
' A; <?H 


tii 

a iH ' 0 rfH 


G 


( 


r)a> 

Ji <9H + 


!c 


*n\ 

i <?H ) ' 


The quantities appearing in the formulas for the perturbations are given in Tables 1-5 and in 
the section on mean motions. 
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Appendix 

Nomenclature 


A - presentation area of a satellite 

A, - coefficients of a Fourier series describing the periodic part of a disturbing function 

A t - coefficients of the Fourier series in the solar radiation pressure disturbing function 

A ffl . - coefficients of the Fourier series in the long period lunar disturbing function depend- 
ent on the moon’s node 

A a (e) - factors of coefficients of the Fourier series which depend on eccentricity of the 
satellite 

A i (i) - factors of coefficients of the Fourier series which depend on the inclination of the 
satellite’s orbit 

A @i - coefficients of Fourier series in the long period solar gravitational disturbing function 

A ci - coefficients of Fourier series in the long period lunar gravitational disturbing function 
independent of the moon’s node 

A Gi - coefficients of Fourier series in the short period solar gravitational disturbing function 
a - semi-major axis of a satellite 

a e - mean value of earth’s radius— (6378.166 km currently in use) 

E - eccentric anomaly of a satellite 
e - eccentricity of a satellite 

e G - eccentricity of sun’s orbit (0.01675184) 
e c - eccentricity of moon's orbit (0.054900489) 

F - radiation pressure constant -4.63 x 10“ 5 dynes per cm 2 
f - true anomaly of satellite 

G - Delaunay variable = faa (l -e 2 ) 

G' - that part of Dalaunay variable G which consists of a constant plus long period terms 
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G* - mean value of G (constant) 

H - Delaunay variable = G cos i 

H' - that part of Delaunay variable H which consists of a constant plus long period terms 
H" - mean value of H (constant) 

i - inclination of satellite orbit plane to earth's equator as a subscript- i indicates a sum- 
mation index 

i ee - inclination of moon's orbit plane to ecliptic ~5.°1 453964 
i £ - inclination of moon's orbit plane to earth's equator 

Jj - earth's second zonal harmonic (1.0822 x 10* 3 ) 

J 3 - earth's third zonal harmonic (-2.285 x 10"*) 

J 4 - earth's zonal harmonic (-2.123 x 10"*) 
j j - integer associated with argument of perigee 
kj - integer associated with right ascension of node 

L - Delaunay veriable = /ila 

L* - that part of L which contains no periodic terms 
L* - equal to L* 
i - mean anomaly of satellite 
i - mean motion of mean anomaly i 

V - that part of mean anomaly consisting of a linear function of time V plus long period 
terms 

l" - that part of l which is a linear function of time 
m - mass of a satellite 

m 0 ' - ratio of mass of sun to combined mass of sun and earth (0.999997) 

m c ' - ratio of mass of moon to combined mass of moon and earth (0.012150668) 

n G - mean motion of sun 0.98560027 degrees per day 
n e - mean motion of moon 13.064999 degrees per day 
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Q(a) - function of semi-major axis a containing also the disturbing force constant 


Q' - disturbing force constant (e.g. for Radiation Pressure Q' = F(A/m )) 

R - disturbing function 

R p - long period part of Radiation Pressure disturbing function 
R g - secular part of a disturbing function 
R 0 ' - short period solar gravitational disturbing function 
R c - lunar gravitational disturbing function 
R 0 - long period solar disturbing function 
r - radius vector from center of earth to satellite 
r ' - radius vector from center of earth to sun 

5 - determining function 
S T - periodic part of s 

- angle between r and r ' 

T - trace of moon’s orbit on celestial sphere 

x , y , z - rectangular coordinates of satellite 
x # , y\ z' - rectangular coordinates of sun 

a. - arguments of trigonometric terms of disturbing function 
a. - mean motion of a. 

a. - arguments appearing in long period radiation pressure disturbing function 
a.' - arguments appearing in short period solar gravitation disturbing function 

e - obliquity of the ecliptic 
\ Q - longitude of the sun 

- k 2 M k = Gaussian constant M = mass of the earth 

6 - argument of perigee of satellite 
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lo - mean motion of w 


w" - linear function of time spearing in o> 
w' = + long period terms 

fi - right ascension of node of satellite 
n" - linear function of time appearing in n 
0' = n u + long period terms 
r - linear combination of «y and 0 
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